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Abstract 

The coherence properties of phase fluctuating Bose-Einstein condensates are studied both theo- 
retically and experimentally. We derive a general expression for the iV-particle correlation function 
of a condensed Bose gas in a highly elongated trapping potential. The second order correlation 
function is analyzed in detail and an interferometric method to directly measure it is discussed 
and experimentally implemented. Using a Bragg diffraction interferometer, we measure intensity 
correlations in the interference pattern generated by two spatially displaced copies of a parent con- 
densate. Our experiment demonstrates how to characterize the second order correlation function 
of a highly elongated condensate and to measure its phase coherence length. 
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I. INTRODUCTION 



Among the various topics related to the exciting field of Bose-Einstein condensation 
(BEC) jl|, the analysis of coherence properties of degenerate Bose gases has attracted major 
interest. Coherence plays a key role in the understanding of the fundamentals of BEC, 
and has a crucial importance for many promising BEC applications, such as matter wave 
interferometry, guided atomic beams, and atom lasers. The coherent character of trapped 
3D condensates well below the BEC transition temperature T c has been confirmed by several 
experiments, using interferometric DO and spectroscopic methods |4(. 

However, recent theoretical and experimental developments have shown that phase coher- 
ence is far from being an obvious property of BEC. In particular, a phase fluctuating BEC _at 
equilibrium has been theoretically predicted in one-dimensional [5(, two-dimensional 



el 3, 



and even in highly elongated, but still three-dimensional |8( trapped Bose gases. Inter- 
estingly, in these cases the density distribution does not differ from the usual BEC profile, 
since density fluctuations are largely suppressed by the repulsive mean-field potential. These 
systems are commonly called quasicondensates. Phase fluctuations can be induced either 
by quantum 3 or by thermal fluctuations For typical experimental temperatures 

quantum phase fluctuations can safely be neglected as long as the system remains in the 
weakly-interacting regime 11]. The amplitude of phase fluctuations, therefore, depends 
strongly on temperature and trapping geometry. In this sense, a nearly phase coherent BEC 
in a highly elongated trap can only be achieved far below T c , imposing severe limitations on 
experiments in constrained geometries. Phase fluctuating BECs have been the subject of 
recent theoretical efforts, including the development of a modified mean-field theory valid 
in all dimensions and all temperatures below the critical point 3, (3 > the analysis of dy- 
namic correlation functions , and the extension of Bogoliubov theory to low-dimensional 
degenerate Bose gases [15^ . 

The phase fluctuating nature of highly elongated BECs was first experimentally demon- 
strated in Ref . 3| • During the ballistic expansion, phase fluctuations transform into density 
modulations. The appearance of phase fluctuations and their statistic nature were studied 
and the dependence of their average value on experimental parameters was characterized 
Moreover, the results obtained from measurements of the energy released during 



the expansion confirmed the absence of density fluctuations in the trapped cloud 3, Q ■ 
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Recently, the physics of quasicondensates has been studied by means of Bragg spectroscopy, 
showing that the existence of phase fluctuations leads to an observable broadening of the 



momentum distribution 



19l l20j . A further experiment has analyzed the phase coherence 



2l|. 



length of non-equilibrium BECs by means of a condensate-focusing technique 

In this paper, we present the theoretical foundation of our studies on coherence properties 
of phase fluctuating condensates. We analyze the behavior of the second order correlation 
function for our experimental conditions and provide a detailed discussion of the experi- 



mental technique used in Ref. [22| to measure it. This technique is based on the analysis 
of the density correlations in the interference pattern generated by a matter wave B ragg 
interferometer. In analogy to the original Hanbury-Brown and Twiss experiment j^l 24 1, 
our method is used to extract the phase coherence length of the degenerate Bose gas from 
density correlation measurements. 

This paper is organized as follows: In Sec. II, we briefly review the theory of phase 
fluctuating Bose-Einstein condensates in 3D elongated traps |8] and analyze the evolution 
of the phase pattern during the ballistic expansion. The knowledge of the free dynamics 
of the phase is important to closely model the BEC evolution during the measurement 
process. In Sec. Ill, we study the coherence properties of the condensate and derive a 
general expression for the iV-particle correlation function of highly elongated 3D BECs. In 
Sec. IV, the experimental technique used to measure the second-order correlation function 
and the phase coherence length of the condensate is reviewed in detail. 

II. PHASE FLUCTUATING CONDENSATES 

n 

In this section, we present the phase operator of a highly elongated condensate [8| and 
develop an analytic description of the ballistic expansion of the fluctuating phase pattern. 
These results, when combined with the free evolution of density modulations presented in 



provide a full understanding of the order parameter dynamics during the time-of- 

fiight. 
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A. Phase operator 



In the following, we consider a cylindrically symmetric condensate in the Thomas-Fermi 
regime, where the repulsive mean-field interaction exceeds the radial (hu p ) and the axial 
(hu z ) trap energies. At T = 0, the density profile has the well-known shape no(p,z) = 
Wom(l — P 2 /R 2 ~ z 2 /L 2 ), where nom = p/g denotes the maximum density of the condensate, 
pi is the chemical potential, g = Anh 2 a/m the interaction constant, m the atomic mass, and 
a > the scattering length. Under the condition u p ^> u z , the radial size of the condensate, 
given by the Thomas-Fermi radius R = is much smaller than the axial size, 

which corresponds to the Thomas-Fermi length L = (2/z/mii;;?) 1 / 2 . 

Due to the repulsive mean-field energy, density fluctuations are strongly suppressed in a 
trapped BEC. Therefore, the field operator describing the condensate can be written in the 



form ip(r) = J n (r) exp(i0(r)), where the phase operator is defined by (see e.g. Ref. 

oo 

0(r) = [4n o (r)]- 1 / 2 5:/;(r)a J + h.c.. (1) 

3=1 

Here aj represents the annihilation operator of the quasiparticle excitation with quantum 
number j and energy 6j] fj" = Uj + Vj is the sum of the excitation wavefunctions Uj and Vj, 
obtained from the corresponding Bogoliubov-de Gennes (BdG) equations. The low-energy 
axial modes, which are responsible for the long wavelength axial phase fluctuations, have 



the energy spectrum e±_ = hu z yjj(j + 3)/4 [2y|. The wavefunctions of these quasiparticle 
modes have the form 



+ (j + 2)(2j + 3)^ (r) fz\ 

4n(j + l)R^Le, j Kl) ' [) 

where Pj are Jacobi polynomials. Equations (JTJ) and (J2J) show that the phase operator 
only depends on the axial coordinate z. In sec. Ill, we analyze the coherence properties of 
the condensate by studying the correlation functions of the operator tp(r). 



B. Evolution of the phase fluctuating pattern 

Starting from the results presented in Refs. 12(, we analyze the evolution of phase 
fluctuations during the free expansion of the degenerate Bose gas. Since the trap is highly 
elongated, we can assume the condensate as an infinite cylinder, and use the local density 



approximation. The time-of-flight dynamics of the order parameter is described by the 



scaling law 



27 



tp{P, z, t) = x ^ e 2BA » e , 



(3) 



where {m\ p /2fi\ p )p 2 is the quadratic phase associated with the expansion dynamics, A 2 (t) = 
1 + uj 2 t 2 is the scaling coefficient, t = J dt' / \ p {t') 2 is the re-scaled time, and p = p/\ p (t) is 
the re-scaled radial coordinate. Let k = tJuq be the solution of the following equation 



_ft_V 2 , ^ P ~2 

2m p 2 



p +g\K \ 



p 



/t = 0. 



(4) 



If we define k = v / raexp(i0), with n = n + 5n, and substitute the scaling law of Eq. (J3J) into 
the corresponding Gross-Pitaevskii equation (GPE), after linearizing in Sn and we obtain: 

d{8n) icp h d 2 



dt 

d(n <t>) 
dt 



A 2 (£) m dz 2 



("00), 



i{Sn/n ) h d 2 

+ ——{on) 



4m cfe 2 



grip 
h\ 2 {t) 



(Sn) 



(5) 
(6) 



where £ = — (7t/ra)[noV~ + VpWoVp]. The first term on th e rig ht hand side of Eq. (JHJ) can 
be neglected in the Thomas- Fermi regime. Following Ref. 26[, we average over the radial 
coordinates. Let nj be the radially-integrated unperturbed density, and 5nj the radially- 
integrated density fluctuations. From Eq. © we obtain: 



[Z,T 



%0) + 



o 2 



:/2 



c 



8A 2 C dz 
5m(z,T') , 

—dr, 



T 5m(z, t') 
o ni(z,T') 



dr' 



(7) 



with r = ujpt, z 



2 Jo A 2 (r') m(z,r') 
z/L, ( = p/hu p , and A = uj p /uj z . Equation (JJJ) can be evaluated from the 



known expression [1( 



5ni(z,r) 



' 2 ) sin 



. -z 2 



nj(z,r) 

where bj = (ej/hup) 2 , dj = bj/( and 

= r (j + 2)(2j + % 1 1/2 (a 3 + a*) 
Cj [AnR 2 Lej{] + 1) J 2 

Q!j and a* are random variables with a zero mean value and (|<x,'| 2 ) 
ing the occupation of the quasiparticle mode j. Near the trap center, Snj/nj 



(8) 



(9) 



Nj, Nj be- 



J2j c jPj\z) sin(ajT)r bj , and hence 



0( 



2,r = 



0(5,0) 




(10) 



For large A and sufficiently short times-of-flight, the significant contribution to the phase 
fluctuations is due to the modes j such that r << A 2 £/[j(j+3)/4], and bj = j(j+3)/4A 2 <C 1. 
Then, using Eq. ((TJ) for 0(5,0), we obtain: 



The second term in the brackets is the correction to the phase contribution of the j-th 
mode due to the ballistic expansion. For typical times-of-flight (tens of milliseconds), this 
correction term is very small (~ 10 -5 ) and the phase pattern can be assumed as completely 
frozen. Using Eq. (fTUJ) . we have verified that, for our typical experimental parameters (see 
Sec. IV-C), the phase change due to the free evolution of the condensate is less than 7r/10. 

III. CORRELATION FUNCTIONS OF A PHASE FLUCTUATING CONDEN- 
SATE 

The coherence properties of a condensate are described by the correlation functions of the 
field operator ip. The importance of correlation functions becomes clear if we consider that 
most experimental signals can be modelled by using this formalism. For example, the first 
and second order correlation functions, describing the single-particle and two-particle corre- 
lation properties of the system, are connected to the visibility of fringes in an interference 
experiment and to the two-body collision rate in the condensate, respectively. 

As discussed in Ref. jjjjjj, the single-particle correlation function of a highly elongated 
degenerate Bose gas can be expressed in terms of the mean square fluctuations of the phase: 




(11) 




(12) 



where 50(r 1; r 2 ) = 4>(ji) — 0(r 2 ) depends directly on the phase operator given in Eqs. (Q). 
At equilibrium, the population of the j-th quasiparticle mode, (a!-Oj), is a random variable 
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with mean value Nj, given by the Bose- Einstein distribution function. The appearance of 
phase fluctuations is a stochastic process governed by the temperature T of the system. 
Since individual realizations are not predictable, we average over an ensemble of identically 
prepared condensates in thermal equilibrium at temperature T. This average is indicated 
by (. . .)t- When k^T ^> hu z is the Boltzmann constant), the population of the j-th 
mode is Nj ~ ksT/ej, and the thermal average of the mean square fluctuations of the phase 
becomes 

<[*fe, *a)] a >T = Sl(T)f 1 \z 1 /L,z 2 /L), (13) 



where 



and 



32MgT 



1 ~ (j + 2)(2j + 3) 



(i.i) PM p(i,i)Pa 



iVo indicating the number of atoms in the condensate fraction. The first order correlation 
function of the degenerate Bose gas is defined by (see e.g. 2^ |) 

W (r r) = (^ t (r 1 )^(r 2 )) T 

' ((^t( ri )^(r 1 ))T^t(r 2 )^(r 2 )) T )V2- 

According to Eqs. (fT2*|) and f(13j) . this results in 

#(^ 2 ) = exp{-5l(T)/( 1 )(z 1 /L,^ 2 /L)/2}. (17) 



For \zi\, \z 2 \ <C L, using the asymptotic expression of the Jacobi Polynomials |3Jj, and 
summing over the different modes in the continuous limit, one obtains an approximated 
formula for the /W function valid around the center of the condensate j^|: 

f^\z 1 /L,z 2 /L) = \z 1 -z 2 \/L. (18) 

In that case, 

g?(z 1 ,z 2 ) =exp{-6 2 L (T)\z 1 -z 2 \/2L}. (19) 
This result suggests the introduction of the phase coherence length of the condensate 

L * = Wry < 20) 

defined as the distance at which the first order correlation function decreases to 1/y/e. The 
approximate formula shown in Eq. (JTSJl can be extended to describe the behavior of the 



function far from the center of the condensate. For 5\{T) ^> 1, the coherence length is 
small compared to the axial size L, and the system is well described by means of the local 
density approximation 16, 17, 2(3]. As pointed out in Ref. j^J, this limit is equivalent to the 
use of the approximate formula for the Jacobi polynomials with large j (30]. Equation (|15|) 
can thus be written in the form 

\z 1 - z 2 \/L 



fU( Zl /L,z 2 /L) 



(21) 



[l-(z 1 + z 2 y/(2L)^- 
generalizing the result obtained in Eq. (|18|). 

We use a similar approach to calculate the two-particle correlation function of the con- 
densate. Introducing the operator 5^(^(ri, r 2 , r 3 , r 4 ) = 4>(ri) + 0(r 2 ) — 4>(r 3 ) — </>(r 4 ), we 
obtain 
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(^(rO^^MrsM^)) = LI VW) exp {-([5 (2) 0( ri , r 2 , r 3 , r 4 )] 2 )/2}. (22) 



i=l 



Using Eq. (0) for the phase operator, a straightforward calculation yields 

ftxWZf M2\ ^ (j + 2)(2j + 3)/i Ar 

([5 { >(j)( Zl , z 2 , z 3 , z A )\ ) = ^ __ T _ r ___JV J 



3 tl 15(i + l)e,-7V 
x 



In the limit k^T ^> Uu z , the thermal average of Eq. (|23|) gives 

([^(zi, za, z 3 , z A )] 2 ) T = 6l(T)fW( Zl /L, z 2 /L, z 3 /L, z A /L), 

where 



(23) 



(24) 



fW( Zl /L,z 2 /L,z 3 /L,z 4 /L) = f^(z 1 /L,z 3 /L) + f^(z 2 /L,z i /L) 

-fU(z 1 /L,z 2 /L)-fW(z 3 /L,z 4 /L) 
+f^( Zl /L lZi /L) + f^\z 2 /L, Z3 /L). 



(25) 



Thus, the two-particle correlation function can be expressed as a product of one-particle 
correlation functions. Equations (fT%|) and (|21j) can be used to derive simplified expressions 
for the function, valid in the limit \zi\ -C L [i = 1,...,4) and in the local density 
approximation. Figure ^ shows the dependence of calculated in 

d + s _ —d — s 



Zl 



22 



2.3 



-d + s _ d — s 

— n ! ^4 = 



(26) 




FIG. 1: f( 2 \zi/L, Z2/ L, Z3/L, z^/L) as a function of s > 0. The complete expression in Eq. 1)25(1 
(solid line) is compared with the approximated formulas derived from Eqs. and (|21j) . valid in 
the condensate center (dotted line) and in the local density approximation (dashed line). The inset 
shows /( 2 ) for different values of d > 0. 

as a function of s > 0. The full expression of can be compared with the two approximated 
formulas, the first valid in the condensate center, the second valid in the local density 
approximation. The inset of Fig. Q shows the same curves for different values of d > 0. This 
choice of variables follows the particular experimental realization. In Sec. IV, we demonstrate 
how these curves can be measured in a matter wave interferometry experiment. There, d 
is the displacement between the two interfering condensate copies, and s is the separation 
between the positions in the interference pattern at which the particle densities are evaluated. 

A qualitative understanding of the behavior shown in Fig. [T] is possible if we consider 
that 

([S^}(z lt z 2 , z 3 , z,)f) T = ([S}(z lt z 3 )] 2 ) T + ([64>(z 2 , z A )f) T 

+2(50(z h z 3 )50(z 2} z 4 )) T . (27) 

The first and the second term are the thermal averages of the operator (<50) 2 calculated in 
(z\, z 3 ) and in (z 2 , £4); the last term is proportional to the correlation function of 5<p at the 
same coordinates. For a fixed displacement d, when the examined positions are close to the 
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condensate center (d, sCi), the first two terms of Eq. (|27|1 do not depend on the separation 
s. However, as s rises from to d, the third term increases from —2([6(j>(zi, ^3)] 2 )r (complete 
anticorrelation) to its maximum value 0, resulting in an uncorrelated phase difference for 
every s > d. In the interval < s < d, the function depends linearly on s with slope 2. 
The second order correlation function is defined as 

#( ri ,r 2 ,r 3) r 4 ) = , ^fr)^^)^^)^^)),. 

((#(r 1 )^(r 1 )} T ...(#(r 4 )^r 4 )} r )i/2 

Substituting Eqs. (J22D and (j2IJ) in Eq. we obtain: 

gf{z^z z ,zti = exp{-5l(T)j< 2 \z 1 /L,z 2 /L,z 3 /L,z,/L)/2}. (29) 
Note that, due to the suppression of density modulations, the normalized density correlation 

(2) 

function of the trapped condensate is constant: (zi, Z2, z 2 , z\) = 1. 

The calculation we have described for the second order correlation function can be ex- 
tended to obtain a general expression for the iV-th order correlation function. Defining the 
operator 

^({r;}^!,...,^) = 0(n) + . . . + - 0(rjv+i) - ... - 0(r 2JV ), (30) 

the iV-particle correlation function is given by 

N 



$t(n) . . . ft(r N )$(r N+1 ) . . . i,(v 2N )) = J] V"o(ri) exp {-([^^({ r! } t=1 ,..., 2JV )] 2 )/2}. (31) 

i=l 

In general, the thermal average of the operator (5^)0) 2 can be written in the form 

<[^)0({r^ =1 ,..., 27V )] 2 > T = 8l(T)f( N \{z i /L} i=1 _ 2N ). (32) 

The function, depending on the Jacobi polynomials P^- 1 ' , can be expressed as a com- 
bination of f^ 1 ' functions: 

/ (J "(WW = E ^' m} / (1) (7, 7O , (33) 

where the coefficient 'P^'™'' is defined as 

,. , , hi if I < N < m 

— 1 if I, m < N or I, m > N 
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The N-ih order correlation function is given by 

(% r |_ _ (j^(ri) . . . fo(T N )4>(T N+1 ) . . . ^(r 2jV )) T 5 

""' (#(r 1 )Vi(r 1 )) T ...(Vit( r2JV )^( r2iV )) T )i/2 

and, from Eqs. (j3*Tj) and (jS^J), 

<7?°({^}i=i,... l2 Ar) = exp{-5l(r)/W({^/L} i=1) ... j2JV )/2}. (36) 

This general result shows that the spatial correlation function of phase fluctuating con- 
densates is completely characterized by the parameter 8\(T) and, therefore, by the phase 
coherence length L&. 



IV. INTERFEROMETRIC MEASUREMENT OF THE SECOND ORDER COR- 
RELATION FUNCTION 

The coherence of a matter wave can be studied by using interferometric methods. How- 
ever, as standard interference experiments measure the first order correlation function of 
the field operator ip, they are very sensitive to phase noise introduced by the experimental 
apparatus. The method presented here is analogous to the original Hanbury-Brown and 
Twiss experiment j3, |3| in which the spatially resolved second order correlation function 
gr^(ri, T2, i"i) of a light source is obtained from intensity correlation measurements. As 
discussed before, for a highly elongated BEC g^\z 1} z 2} z 2 , zi) = 1. This result suggests 
that a simple measurement of density correlations in the condensate is not sufficient to 
describe the coherence properties of the sample. Nevertheless, by measuring density corre- 
lations in the interference pattern generated by two spatially displaced copies of a parent 
BEC, it is possible to correlate the field operator ip at four different positions and extract 
gT (zi, Z2, Z3, Z4J directly. Compared to standard interference experiments, the main advan- 
tage of this technique is the intrinsic stability of the density correlation measurement against 
variations of the global phase between the interfering condensates. 

In this section, we show how a matter wave Bragg interferometer can be used to character- 
ize the second order correlation function of the condensate and measure its phase coherence 
length. 



11 



A. Interferometric scheme 



Our interferometric sequence is shown in Fig. |21 The condensate is released from the 
magnetic trap and expands freely for 2 ms. This short time-of-flight is important to lower 
the density, thus reducing s-wave scattering processes occurring during the Bragg diffraction 



of the condensate |3l(. The interrogation sequence consists of two it/2 Bragg pulses. Each 
pulse is composed of two counterpropagating laser beams of wave number k, detuned from 
the atomic transition. The first Bragg pulse splits the condensate in the two momentum 
eigenstates \2hk) and |0) along the axial direction (z). After a time At, a second 7r/2 pulse 
splits the condensates again, creating two interfering copies in each momentum state. The 
time interval At between the two pulses sets the spatial overlap, d = 2hkAt/m, between 
the interfering BECs at the output ports of the interferometer. The relative phase of the 
two counterpropagating Bragg beams is externally controlled by an electro-optic modulator 
(EOM) and can be changed between the two pulses. This allows us to imprint an extra 
phase (p which can be precisely tuned. Control of the EOM phase is crucial for our method, 
as described in Sec. IV-B. 

Using the results derived in Sec. II, the atoms detected in the output port A (Fig. |2J), 
after a total time-of-flight t, are described by the order parameter 

V(r, d, t) = -y/v&t) + ^y/ri(T,t)exp{i[5<t>(z, z' , t) + a(z, z' , t) + /3{z, z') + 7 (d)]}, (37) 

where r' = r — dz and r)(r,t) is the time-evolved density profile normalized to the total 
number of atoms in the parent condensate. The relative phase between the interfering 
condensates contains several contributions. 54>(z,z',t) = 4>(z,t) —<f)(z',t) describes the 
phase difference between z and z' that evolves from the phase fluctuations in the parent 
condensate. The term 

a(z,z',t) = ^-y-z' 2 ) (38) 

represents the non-uniform spatial phase profile developed during the mean-field-driven ex- 
pansion. The mean-field gradient between the interfering BECs is responsible for a force 
repelling the centers of mass of the two clouds. This effect is described by the phase term 

PM = ^(z + *), (39) 

proportional to the relative repulsion velocity 8v between the interfering condensates j^. 
After the first Bragg pulse the relative phase of the atoms in the \2hk) momentum state 
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FIG. 2: a) Matter wave Bragg interferometer. The condensate is released from the magnetic trap 
and evolves freely for 2 ms. The sample is interrogated by the first ir/2 Bragg pulse which splits 
the parent BEC in two copies with momenta and 2hk. After a time At, the second ir/2 Bragg 
pulse splits the condensates again and allows them to interfere. The time interval At defines 
the displacement d between the two interfering condensates, b) A typical line density profile 
at the output ports of the interferometer. The distance between the two autocorrelated copies 
(d = 46 /urn) is comparable to the phase coherence length of the parent condensate (L^ = 43 /xm).In 
the schematic of the matter wave Bragg interferometer, the distance d has been exaggerated for 
clarity. 
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evolves with a characteristic frequency < 5p r ggg, given by the detuning of the lasers from the 



resonance of the two-photon transition 



33[. Therefore, the last term, 



TYld 

7(d) = S Biagg At + ip = 5 Bragg — + <p, (40) 

represents a global phase depending on the detuning from the Bragg transition and the 
externally controlled phase ip. 

The density of atoms at the output port A of the interferometer is given by 

I(v,d,t)= ± v (r,t) + ^r l (r',t) 

+ ^ V (r,t) V (r',t) cos[5<P{z, z', t) + a(z, z' , t) + (3{z, z') + 7 (d)]. (41) 

The presence of strong phase fluctuations alters the interference pattern generated by the 
two autocorrelated condensates. In fact, when d ~ the phase term 5(f) can be comparable 
to 7T, modifying drastically and in an unpredictable way the position and the spacings of the 
interference fringes. 



B. Method 

Starting from Eq. (}4~Tj) . we want to calculate the density correlation function of the inter- 
ference pattern, for an ensemble of identically prepared condensates at a given temperature 
T, averaged over all the global phase values (p. This averaging process is indicated by the 
symbol (. . )t,^- It is therefore important that the phase delay (p induced by the EOM is uni- 
formly changed between and 2ir. In Sec. II, we have shown that, for typical times-of- flight 
(tens of milliseconds), the evolution of the fluctuating phase of the condensate is basically 
frozen. This allows us to neglect the time-dependence of 5<f)(z, z' ,t). We also neglect the 
contribution of density modulations induced by the initial phase pattern on the Thomas- 
Fermi profile of the condensate. The validity of this approximation is verified below. Under 
these assumptions, we calculate the normalized density correlation function 

7 ( 2 )( ri ,r 2 ,d,t) = ((h - (h) T , v ){l2 - {I2)t, v ))t, v (42) 

where = I(r lj2 , d, t). After a lengthy but straightforward calculation, the averaging 
process gives 

^ 2 \zi,z 2 ,d, t) = cos 




d + 5v) {z x — z 2 ) 
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x exp[-6 2 L (T)f^( Zl /L, (z 2 - d)/L, ( Zl - d)/L, z 2 /L)/2}. (43) 

7^ 2 ) (zi, z 2 , d, t) results from the product of two different terms: the first is a periodic function, 
whose argument is the contribution of the mean-field energy to the phase profile (ballistic 
expansion and relative repulsion between the interfering condensates); the second is an 

(2) 

exponential term which corresponds to the g T function of the parent phase fluctuating 
condensate. The decay constant of this function is given by the phase coherence length of 
the condensate (see Eq. I2U|) . 

From the experimental point of view, the averaging process described above is equivalent 
to the following procedure: The radially integrated density profile / = I(z, d, t) at the output 
port A of the interferometer is measured for different values of the global phase, uniformly 
distributed in the range < <p < 2ti; then the average value {I)t,<p is calculated and used to 
determine I — {I)t,<p for each experimental realization. These profiles, averaged according 
to Eq. give a measurement of ^ 2 \zi, z 2 , d, t). We evaluate the density correlations 

as a function of the separation s — z 2 — z\. For simplicity, we choose symmetric positions 
around the center (z = d/2) of the interference pattern in the output port A. The positions 
in Eq. are defined such that Z\ = z\ and z 2 = z 4 (see Fig. EJ). The method described 
here allows us to characterize the dependence of the correlation function 

m ( \ z 



n\x d + 5v ) s 



7 ( - 2 )(s, d, t) = cos 

x exp[-6 2 L (T)f^\z 1 /L,z 2 /L,z 3 /L,z 4 /L)/2] (44) 



on the separation s for any fixed displacement d between the interfering condensates. 



C. Experimental results and numerical simulations 

We perform the experiment with 87 Rb condensates in the F = 1, mp = — 1 state. The 
atoms are confined in a highly elongated magnetic trap with cylindrical symmetry, the long 
axis lying in the horizontal plane. The confining potential has an axial frequency uj z = 
2tt x 3.4 Hz and a radial frequency u p which is varied between 2ir x 300 Hz and 2ix x 380 Hz. 
Further details on the experimental apparatus can be found in 18]. After the BEC formation, 
we let the system thermalize in the magnetic trap for typically 4 s in presence of radio 
frequency shielding t 34] . That time is important to reach an equilibrium condition in which 
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any quadrupole oscillation has been damped down. As shown in Fig. El our matter wave 
interferometer consists of two ir/2 Bragg diffraction pulses. Each of them is composed of two 
counterpropagating laser beams, detuned by about 3 GHz from the atomic transition. This 
detuning suppresses spontaneous scattering of photons during the interrogation time. The 
Bragg pulse duration of 100 /is is sufficiently short not to resolve the momentum distribution 
of the atoms in the condensate and long enough to avoid higher order Bragg diffraction 
processes. A fixed frequency difference is set between the two counterpropagating beams to 
match the Bragg condition. The condensate is released from the magnetic trap and after 
2 ms of time-of-flight is probed by the two-pulse sequence of the interferometer. The atomic 
cloud is detected after the ballistic expansion by resonant absorption imaging. 

Figure |2b shows a typical line density profile of an interference pattern where the distance 
between the two autocorrelated copies (d = 46 /mi) is comparable to the phase coherence 
length of the parent condensate (L^ = 43 /mi). Because of the stochastic nature of phase 
fluctuations, the fringe spacing is not regular and differs in each experimental realization. 
This experimentally demonstrates that the fluctuating phase of the condensate can signifi- 
cantly change on distances comparable with the phase coherence length of the sample. Even 
if each single image shows high contrast, the interference pattern is completely washed out 
when we average a significant number of realizations. 

The results of standard interference experiments are related to the correlations of the 
wavefunction and therefore are very sensitive to phase instabilities. Figure El shows the 
interference signal obtained by measuring the number of atoms in an interval of width 
0.2 x L around the center of the interference pattern (z = d/2) at the output port A, 
as a function of the global phase ip controlled by the EOM. This signal is normalized to 
the corresponding number of atoms in the parent condensate. The two plots correspond 
to different displacements d between the interfering condensates. A small displacement is 
related to a short time interval between the two interrogation Bragg pulses. In that case, 
the contribution of phase fluctuations and the effect of technical phase noise introduced by 
the experimental apparatus are both negligible. Therefore, according to Eq. (|41jl . when 
d <C Lfi and At is small compared to the characteristic time stability of our Bragg pulses, 
the normalized signal oscillates sinusoidally with high contrast. For d approaching L^, the 
random phase introduced by the phase fluctuations washes out the oscillation. If external 
disturbances can be neglected, the contrast of the oscillations is directly related to the 
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FIG. 3: The number of atoms measured in the interval d/2 — 0.1 x L < z < d/2 + 0.1 x L, around the 
center of the interference pattern detected at the output port A, is plotted as a function of the phase 
(p controlled by the electro-optic modulator. The signal is normalized to the corresponding number 
of atoms in the parent condensate. The two sets of data correspond to different displacements d 
between the overlapping condensates. The solid line is obtained by fitting the experimental data 
with a sinusoidal function. The measurements refer to condensates with about 3 x 10 5 atoms, a 
typical axial size of L = 180 /im and a temperature T = 170 nK . 



first order correlation 
external disturbances 



unction g^ 1 ' at a given displacement d. However, as d increases, the 
35^ also increase and produce a random phase noise which destroys 
the oscillating behavior and hides the effect of phase fluctuations on the detected signal. 

This problem can be solved by using the method described in Sec. IV-B. The measurement 
of intensity correlations, in combination with the subsequent averaging process, has the 
major advantage of being insensitive to technical phase noise introduced by the experimental 
apparatus. Figure 0] shows the correlation function 7^ 2 - ) (s, d, t) extracted from a set of 29 line 
density profiles corresponding to 5.0 x 10 5 condensed atoms at a temperature T = 216 nK, 
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FIG. 4: Circles: Correlation function 7 2 (s, d, t) extracted from a set of 29 line density profiles. The 
data correspond to samples with 5.0 x 10 5 condensed atoms at a temperature T = 216 nK, detected 
after a total time-of-flight t = 37 ms. The displacement between the interfering BECs is d = 35 p, 
The bars on the experimental points represent the statistical errors. Crosses: Numerical simulation 
which takes into account the time dependence of the fluctuating phase and of density modulations, 
modelled on the experimental parameters. Solid line: Fit to the experimental data using the model 
function of Eq. (|45[). Dashed line: Second order correlation function g^, (s,d) = g?p (21,^2,^3,^4) 
extracted from the fit to the experimental data. The phase coherence length of the sample is 
graphically indicated on the plot. 

detected after a total time-of-flight t = 37 ms. The displacement between the interfering 
BECs is d = 35 fim. The experimental data is compared with a numerical simulation which 
produces random phase patterns according to the experimental conditions and uses Eq. ()37|) 
to describe the evolution of the order parameter. The numerically calculated points shown 
in Fig. |3] are obtained by following the same averaging procedure we have applied to the 
experimental data. This kind of analysis includes the time dependence of the fluctuating 
phase and of the density modulations induced by the initial phase pattern. The solid line 
is the result of a fit to the experimental data. According to Eqs. ()44|) and ()26|). the model 
function 

cos(a • s) exp[-b ■ f®(z 1 /L,z 2 /L,z 3 /L,zt/L)/2] (45) 
contains only two free parameters. The curves clearly show the damped oscillating behav- 
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FIG. 5: Direct comparison between the measured phase coherence lengths and the theoretical 
values, calculated according to Eq. ((2f)|) by using the measured numbers of atoms, temperatures 
and trapping frequencies. The dotted line with slope 1 is used to compare experiment and theory. 
The bars on the plotted points indicate statistical errors. The relative systematic uncertainties on 
the calculated and measured pha se coherence length are 26% and 15%, respectively. This figure 
has previously been shown in [22j. 

ior. The oscillation frequency strictly defines the parameter a, while the damping coefficient 
gives a measurement of 8\{T). From the fitting function it is possible to extract the spa- 
tial dependence of the second order correlation function g^\s,d) = (zi,Z2,z^,Z4) (see 
Eqs. ()29|) and 1)440). The fit on the experimental data gives a phase coherence length 
L^ xp = (57 ± 10) fim, compatible with the expected value L^ 1 = (58 ± 2) /im. The good 
agreement between the experimental data, the numerical simulation and the model function 
of Eq. ()3H|) demonstrates that the free evolution of density modulations and of the fluctuat- 
ing phase pattern do not influence the measurement of the second order correlation function. 
This result justifies the use of Eq. (}4"3|) to model the experimental data and to extract the 
coherence properties of the condensate. In Fig. we show a direct comparison between 
the measured phase coherence lengths in the center of the BEC and the theoretical values 
calculated according to Eq. (}2()j) . by using the measured numbers of atoms, temperatures 
and trapping frequencies. The bars indicate the statistical errors both on the measured 
values and on the theoretical predictions. The dotted line with slope 1 highlights the good 
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quantitative agreement between experiments and theory. 
V. CONCLUSION 

In this paper, we have studied the coherence properties of phase fluctuating Bose-Einstein 
condensates. In highly elongated BECs the thermal excitation of quasiparticle modes can 
significantly reduce the coherence length of the system. Starting from the results of Petrov 
et al. jsj], we have derived a general formula for the A-particle correlation function. The 
second order correlation function has been studied in detail and its limits both around the 
center of the condensate and in the local density approximation have been analyzed. In 
particular, we have discussed a method to directly characterize the second order correla- 
tion properties of the system. An analytic theory that describes the free evolution of the 
condensate phase has been developed to closely model the measurement process. Using a 
Bragg diffraction interferometer, we have measured the density correlations of the inter- 
ference pattern generated by two spatially displaced copies of a parent BEC. This kind of 
measurement allows to correlate the field operator ip of the parent condensate in four differ- 
ent z positions. The averaging process directly gives the second order correlation function. 
The experiment confirms our theoretical predictions and demonstrates a method to measure 
the phase coherence length of the condensate. Compared to usual interference experiments 
this technique has the advantage of being insensitive to the global phase noise introduced by 
the experimental apparatus. The method presented here is in direct analogy to the original 
Hanbury-Brown and Twiss experiment and demonstrates the possibility of using density 
correlation measurements to study the coherence properties of Bose-Einstein condensates. 
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